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ABSTRACT 

Context. Near-Earth asteroid-comet transition object 107P/ (4015) Wilson-Harrington is a possible target of the joint European Space Agency 
(ESA) and Japanese Aerospace Exploration Agency (JAXA) Marco Polo sample return mission. Physical studies of this object are relevant to this 
mission, and also to understanding its asteroidal or cometary nature. 

Aims. Our aim is to obtain significant new constraints on the surface thermal properties of this object. 

Methods. We present mid-infrared photometry in two filters (16 and 22 pm) obtained with NASA's Spitzer Space Telescope on February 12, 
2007, and results from the application of the Near Earth Asteroid Thermal Model (NEATM). We obtained high S/N in two mid-IR bands allowing 
accurate measurements of its thermal emission. 

Results. We obtain a well constrained beaming parameter (rj = 1.39± 0.26) and obtain a diameter and geometric albedo of D = 3.46 ± 0.32 km, 
and pv = 0.059 ± 0.01 1. We also obtain similar results when we apply this best-fitting thermal model to single-band mid-IR photometry reported 
by Campins et al. (1995), Kraemer et al. (2005) and Reach et al. (2007). 

Conclusions. The albedo of 4015 Wilson-Harrington is low, consistent with those of comet nuclei and primitive C-, P-, D-type asteorids. We 
establish a rough lower limit for the thermal inertia of W-H of 60 JnT^ s^°'^ K^^ when it is at r=lAU, which is slightly over the limit of 30 
Jm^^s^^ ^ K^' derived by Groussin et al. (2009) for the thermal inertia of the nucleus of comet 22P/Kopff. 
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1. Introduction 

The asteroid-comet transition object 107P/ (4015) Wilson- 
Harrington (hereafter W-H) is a possible target of the joint 
European Space Agency (ESA) and Japanese Aerospace 
Exploration Agency (JAXA) Marco Polo sample return mission. 
This mission seeks to understand the origin and nature of volatile 
and organic material in the early Solar System and ideally would 
do this by returning samples of a primitive asteroid - i.e. a C, P 
or D asteroid in the Tholen classification system (Tholen and 
Barucci 1 989 1). C, P and D objects are abundant in the middle to 
outer Main Belt; and are believed to be primitive remnants from 
the early Solar System. 

W-H was discovered with a cometary appearance in 1949 
(Cunningham 1950, Fernandez et al. 119971 ). Recovered in 1979, 
and observed at every apparition since, it has always appeared 
as a point source. W-H may have a cometary origin, but its 
cometary nature is not firmly established. Its orbit has an un- 
likely cometary origin, its Tisserand parameter is 3.084 and it has 
only a 4% chance of coming from the Jupiter-family comet pop- 
ulation (Bottke et al. |2002| l. However, the same argument against 
a cometary origin could be made for the active comet 2P/Encke, 
which could have achieved its orbit either due to the influence of 
nongravitational forces (Fernandez et al. 2002) or gravitational 
perturbations (Levison et al. 2006). Both photometric and spec- 
troscopic data suggest classification of W-H within the C com- 
plex (Tholen and Barucci [T989l Chamberhn et aI. |T996] l. 
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Licandro et al. ( |2007b| l suggested a link between W-H, 
(3200) Phaethon (another similar asteroid-comet transition ob- 
ject in the NEA population), and the objects with asteroidal or- 
bits in the main belt that show cometary-like activity (the "Main 
Belt Comets", MBCs, Hsieh & Jewitt 2006), as their colors 
and/or spectra correspond to those of Cb- or B-type asteroids. 
Licandro et al. show that Phaethon's surface composition sug- 
gests it is unlikely that this object has a cometary origin. 

If the mechanism that activates the MBCs is driven by the 
sublimation of ices, that means some asteroids in the outer main- 
belt could have retained ices below their surface and, under cer- 
tain conditions, become "activated asteroids". This supports the 
scenario that the snow line fell somewhere in the asteroid belt 
around 3AU. This suggest that asteroids could have contributed 
to the formation of the oceans on Earth due to impacts on our 
planet (e.g. Drake and Campins 2006 Motd et al. |20071 l. While 
W-H has shown cometary activity, the existence of MBCs shows 
that this does not preclude it from having an asteroidal rather 
than cometary dynamical origin. 

Mid-infrared observations are needed in order to derive an 
accurate size and albedo using the near-Earth asteroid thermal 
model (NEATM). The surface temperature can be used to con- 
strain the thermal inertia . Information about the surface temper- 
ature can be obtained from multi-wavelength filter photometry. 
We present new high S/N mid-infrared broadband photometry 
of W-H which places significant new constraints on its surface 
properties, and we compare these results with previously pub- 
lished mid-infrared data. 
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Date(UT) r (AU) p (AU) q-(°) Fjs.siinJy) Fji.ijmJy) 
QliniQl 4.078 3.520 12.7 1.26±0.08 2.20±0.04 

Table 1. Observing geometry and measured fluxes . 



2. Observations and data reduction 

W-H was observed as part of the SEPPCoN program (Survey 
of Ensemble Physical Properties of Cometary Nuclei, see 
Fernandez et al., I2008I I. with the Spitzer Space Telescope on 
February 12, 2007. Observations started at 07:36:55 UT and 
lasted for 17.14 minutes. We used the peak-up imaging mode 
of the InfraRed Spectrograph (IRS) and obtained images in both 
the "blue" and "red" passbands (monochromatic wavelengths of 
15.769 and 22.327 //m, respectively). The total integration time 
was 94 s for the "blue" and 147 s for the "red" images. The data 
were processed and photometrically calibrated with the Spitzer 
Science Center's IRS pipeline (version 16.1.0). MOPEX soft- 
ware (version 16.2.8) (Makovoz and Khan. l2005T l was used to 
clean, register, and co-add the data. Each exposure was aligned 
onto a common coordinate system, taking into account the mo- 
tion of W-H, and any constant offsets in the background were re- 
moved. The exposures were averaged together, filtering out bad 
pixels, which produced two final images: one each for the blue 
and red peak-up arrays. Our data reduction procedure is identi- 
cal to that followed in our previously described work (Groussin 
et al. I2009I I. The object appeared as a point-source in both im- 
ages. The photometric analysis was then performed by means of 
point spread function fitting . Color correction to the photometry 
was applied based on W-H's color temperature. The observing 
geometry and measured fluxes are in Table [T] where r and p are 
the heliocentric and Spitzer distance respectively, a is the phase 
angle, F15 « and F22.3 are the fluxes measured in the "blue" and 
"red" bands. 



3. Thermal modeling 

The flux of asteroids in this spectral range is dominated by ther- 
mal emission. The measured spectral energy distribution (SED) 
depends on the object's size, composition, and temperature dis- 
tribution. This last term is dependent on several factors, in- 
cluding distance from the Sun, albedo, thermal inertia, surface 
roughness, shape, rotation rate, and spin-pole orientation. Many 
of these are unknown for W-H, therefore fitting our photome- 
try with a thermophysical model (Muller 120071 and references 
therein) is unwarranted. We chose to use the semi-empirical 
Near-Earth Asteroid Thermal Model (NEATM, Harris 1998) to 
fit our 16 and 22 /im photometric data points in order to derive 
the size and albedo of the object. We also use the NEATM to 
derive information on the thermal intertia of W-H. We use the 
NEATM to solve simultaneously for the beaming parameter (77) 
and the diameter {D). The beaming parameter was originally in- 
troduced in the STM to allow the model temperature distribution 
to fit the observed enhancement of thermal emission at small 
solar phase angles due to surface roughness. In the NEATM, rj 
can be thought of as a modeling parameter that allows a first- 
order correction for any effect that influences the observed sur- 
face temperature distribution (such as roughness, thermal iner- 
tia, and rotation). At small phase angles (a <30°) values of 77 < 
1 indicate that the presence of surface roughness enhances the 
temperature with respect to that of a smooth surface, whereas 
77 > 1 indicates substantial thermal inertia. The values of 77 , D 
and the geometric albedo {pv) derived for 4015 W-H using the 
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Table 2. NEATM Best Fit Model of the dat presented in this 
paper, and the results from Kxaemer et al. (2005). 
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Fig. 1. Thermal model fit (solid line) using NEATM to the ob- 
served flux (circles) of 4015 Wilson-Harrington. 

NEATM (see Fig.[TJ are presented in Table |2l Uncertainties in 
the diameter estimated with thermal modeling usually exceed the 
formal errors and are typically 10-15% (e.g., Delbo et al. l2003l l. 
We assumed a thermal infrared emissivity (e) of 0.9. To estimate 
Pv we adopted Hv = 15.99 ± 0.10 and G = 0.15 (MPC 17270), 
where Hy is the visible absolute magnitude and G is the slope 
parameter. The listed error in the geometric albedo includes the 
uncertainty in Hy and the diameter uncertainty propagated into 
the albedo. The error we assume in Hy (0. 1 mag) is due to the 
rotational lightcurve amplitude of 0.2 magnitudes (Osip et al. 
[T995T l. 

4. Discussion and conclusions 

Our results can be compared with previous results. Mid-infrared 
observations of W-H were reported by Campins et al. ( |1995l l. 
Kraemer et al. (2005), and Reach et al. ( |2007i i. 

Campins et al. presented simultaneous J, H, K, and N band 
broad-band photometry and derived two different values of size 
and albedo, a diameter of 2.60±0. 16 km and a geometric albedo 
in the J band of 0.10±0.02 using the STM, and a diameter 
of 3.90±0.25 km and albedo 0.05±0.01 using the Isothermal 
Latitude Model (ILM). Kraemer et al. presented broadband pho- 
tometry of W-H in two filters centered at 8.3 and 14.7 jum, us- 
ing the SPIRIT III instrument on board of the Midcourse Space 
Experiment, and using NEATM derived D, pi/ and rj (their re- 
sults are presented in the second line of Table |2|. Reach et al. 
presented Spitzer Space Telescope images using the 24 /zm fil- 
ter, and reported the measured emission flux in this band, but 
they did not derive either the size or the albedo. 

We first compare our results with those of Kraemer et 
al. since they obtained photometry in two bands, and applied 
NEATM to those data (see Table While their D and pv are 
consistent with ours, their rj value is smaller. Their observations 
were done at a much higher phase angle, so this doesnot fit well 
with the trend of increasing rj with increasing solar phase an- 
gle that has been observed for about 40 near-Earth asteroids so 
far (e.g., Delbo |2004l Campins et al. |2009l l. However according 
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Table 3. NEATM Best Fit Models to previously published data assuming a fixed value of 77 = 1 .39 obtained in this paper, and this r\ 
+ 1-cr values in order to derive a possible range. 



to Kraemer (personal communication) they did not use the 14.7 
ytva flux to derive their reported D and as they concluded that 
this flux was poorly constrained. So the D and py reported by 
Kraemer et al. were obtained by fitting NEATM to the 8.3 yum 
flux alone while using a fixed value of 77. This means that our 
computed value of 77 represents the first measurement of W-H's 
beaming parameter. 

Finally, we used the NEATM to fit only the Kraemer et al. 
flux at 8.3 yim using our computed value of 77 = 1.39 ± 0.26. To 
incorporate the uncertainty in 77, we used 3 different values 77 - 
1.39, 1.13, and 1.65 (= 1.39 ±lcr) respectively. We did the same 
with the Campins et al. and Reach et al. observations. Notice 
that in those cases the phase angle is very similar. In the case of 
Reach et al., we re-analyzed the images and computed the flux, 
and after applying the color correction we obtained a value of 
3.8 mJy at 23.68 //m. Results are shown in Table[3] where CWL 
is the central wavelength of the reported fluxes (Flux) . 

From Table [3] we can conclude that the D, pv and 77 val- 
ues we determine with the Spitzer data are robust and consistent 
with previously published observations. In particular, the albedo 
obtained is consistent with that of cometary nuclei (Fernandez et 
al. |2005 ) and of primitive asteroids. It is interesting to note that 
the 77 value we obtained (1.39) is high at a rather small phase an- 
gle. Notice also that, in all cases, to have a cometary like albedo 
(pv < 0.075 according to Fernandez et al. 120051 1, W-H's 77 should 
be larger than 1.1. 

The analysis of the SEPCoN sample by Fernandez et al. 
(2009 in prep.) report that they have 21 cometary nuclei (in- 
cluding W-H) for which the measured beaming parameter 77 is 
probably robust. While the 1-cr error bar on each value of eta 
is comparable to what we report here for W-H (±0.2 to 0.3), the 
distribution of 77 values is wide, with two-thirds of the group near 
1 .0 and the other third straddling 1 .4. So, the 77 value of W-H fits 
well in the second group of comets and cannot also be used to 
discriminate between the asteroidal or cometary nature of this 
object. 

On the other hand, the value of 77 can be used to estimate the 
value of the thermal inertia from a relation existing between 77 
and the thermal parameter = F /(ecrT^^), where F is 

the value of the thermal inertia, P is the rotation period, e the 
infrared emissivity, cr the Stefan Boltzmann constant, and Tss 
the subsolar temperature. Tss can be derived from the equation: 
(l-A)SQr^^ = TjecrT^^ , where A is the bolometric Bond albedo, 
and 5o is the solar constant at 1 AU. The latter can be calculated 
from the value of pv in Table |2] using the formula A = pv x 
(0.29 -I- 0.684G) then A=0.23 (Bowell et al., [T989] l and assuming 
a slope parameter G=0.15. 

The derived value of 77 of about 1 .4 implies a value of be- 
tween 1.5 and 3.5 according to Spencer (.1990; Fig 5), depending 
on the degree of surface roughness (which is unknown for W-H) 
and assuming observation at a phase angle cc = 0° (a valid com- 
parison considering that our observations were done at a very 
low phase angle, and thus the possible variation of 77 with a is 
much lower than the uncertainties) . This corresponds to a range 



of thermal inertia between 20 and 70 Jm "^A" ' for our ob- 
servations obtained at r=4.078 AU. 

The thermal conductivity in the regolith of an airless body is 
temperature dependent (Keihm, 1984), as is the thermal inertia. 
This temperature dependence of F must be taken into account in 
order to compare our results with those obtained for NEAs (e.g. 
Delbo et al.. i2007l i that are usually observed closer to the Sun (at 
about 1 AU). Under the assumption that heat is transported in the 
regolith mainly by radiative conduction between grains, the ther- 
mal conductivity is proportional to T^^ (Kiihrt and Giese, [19891 

Jakoskv, [l986l ). In this case F oc T^^^ and thus F oc p/'^. On the 
basis of this dependence we find that the range of thermal inertia 
is about 60 and 200 Jlr^'^s^^^^^K^^ at r-\ AU. These values are 
between the typical thermal inertias of km-sized NEAs (Delbo 
et al., 120071 ) and those of cometary nuclei (Groussin et al. |2009l l. 
We caution here that the values above should be taken as indica- 
tive because they assume a spin-pole orientation perpendicular to 
the orbital plane and a -Q° . Unfortunately the spin-pole of W-H 
is unknown, so larger values of F cannot be ruled out. Moreover, 
< 5% lower values of F would have been derived if the morning 
side (at a - 12.7°) of the body was actually observed. 

In conclusion, the albedo of 4015 Wilson-Harrington is low 
(pv-Q-Q59 + 0.011), consistent with those of comet nuclei and 
primitive C-, P-, D-type asteorids. The value of its beaming pa- 
rameter (77=1.39+0.26) is compatible with about 1/3 of the 21 
comet nuclei with accurate 77 values determined by Fernandez et 
al (2009, in prep.). From 77, we derived a lower limit for the W-H 
thermal inertia of W-H of 60 ym^^i '^^/T which is slightly 
over the limit of 30 Jm^^ s^'^*^^ K^^ derived by Groussin et al. 
( 120091 ) for the thermal inertia of the nucleus of comet 22P/Kopff. 
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